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We investigate the nature of the three-mode interaction inside an optomechanically-active micro- 
toroid with a sizeable x^'^^ coefficient. Experimental techniques are quickly advancing to the point 
where structures with the necessary properties can be made, and we argue that these provide a 
natural setting in which to observe rich dynamics leading, for instance, to genuine tripartite steady- 
state entanglement. We also show that this approach lends itself to a full characterisation of the 
three-mode state of the system. 



Over the last several years, the field of optomechan- 
ics has witnessed remarkable progress in experimental 
achievements [1, 2], the chief driving factor behind which 
was the quest to achieve ground-state cooling of a me- 
chanical oscillator. This was achieved, first in by means 
of cryogenics [3], then by electromagnetic means in an 
electromechanical system [4], and finally in an optome- 
chanical setting [5] . Reaching the ground state is a means 
to an end, for it is only when a mechanical oscillator is 
close to, or at, the ground state that its true quantum 
nature shows up unambiguously. This was demonstrated 
clearly in a recent experiment [6] that showed the im- 
balance between the red- and blue-mechanical sidebands 
in the spectrum, a clear signature that the mechanical 
oscillator is behaving in a nonclassical way. 

Turning away from fundamental physics, one would 
like to use quantized mechanical resonators as a resource] 
typical oscillators have decay rates in the sub-kHz 
domain, meaning that the decoherence time 



can be made large compared to the other timescales 
of the system by using cryogenic methods to decrease 




FIG. 1. Schematic diagram of the system, (a) Side view of 
the toroid. (b) Top view. We show the coupling of the cavity 
fields to the field in the waveguide. 
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^th = ^B^env/(^^m), the average number of phonons 
at an environmental temperature Tenv and at mechani- 
cal frequency uo^ (/cb is Boltzmann's constant), as much 
as possible. Clearly, mechanical oscillators with a large 
mechanical frequency, say cjm ^ 27r x 1 MHz, and large 
mechanical quality factor Qm = ^m/ are at an ad- 
vantage in this respect. At the same time, one would like 
the mechanical oscillator to interact strongly with an op- 
tical resonator that has a similarly large optical Q. It is in 
this context that optomechanical toroidal structures [7] 
appear as ideal optomechanical systems. From a techno- 
logical point of view, toroidal structures are also ideal in 
that they minimize the number of moving parts — there 
are no moving mirrors to align — and can be manufac- 
tured monolithically on CMOS-compatible substrates [8], 
pointing the way towards a possible integration with con- 
ventional (opto)electronics in the future. 
The highest-quality optical modes in toroidal structures 
are of the 'whispering gallery mode' type, with the me- 
chanics of total internal reflection ensuring that losses 
are minimized. A recent group of experiments [8, 9] has 
recognized this feature as enabling another technology: 
second-harmonic generation (SHG). Indeed, it turns out 
that toroidal and ring-resonator structures facilitate SHG 
because the phase matching conditions that are necessary 
in any nonlinear optics experiment can be met automat- 
ically by choosing the right doublet of optical modes [9] . 
It is the purpose of this paper to combine these two ideas. 
We shall look at the emergence of nonclassical steady- 
states, e.g., genuinely tripartite entangled states, in the 
three-mode system formed by the two optical modes — the 
fundamental and the second harmonic — and the mechan- 
ical oscillator. Every pair of these three modes interacts 
directly, and we shall see that this results in a competi- 
tion between the two purely optomechanical interactions 
and the second-harmonic generation process. 

This paper is structured as follows. Over the next sec- 
tion we shall introduce the full model Hamiltonian, and 
then proceed to obtain the equations of motion. The 
usual procedure is used to linearize the dynamics, where- 
upon we can concentrate exclusively on Gaussian states 
and present some numerical results. The next section 
discusses state detection using homodyning techniques. 
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after which we conclude our investigation. 



I. PROPOSED MODEL 

The model Hamiltonian we use is a combination of 
the usual optomechanical Hamiltonian and the Hamilto- 
nian description of the SHG process, and describes the 
system shown schematically in Fig. 1. We shall label the 
annihilation operators of the two optical fields ap (funda- 
mental 'F', of frequency uJc) and as (second harmonic 'S', 
frequency 2co'c). These two modes are coupled to a contin- 
uum of modes, represented by the operators a^^ through 
decay rates hz^^s (in practice, /^p ^ i^s [9]), as well as 
to each other through a second-harmonic interaction fre- 
quency X- The mechanical oscillator 'M' is represented 
through its dimensionless quadratures x and p, and is 
characterized by the mechanical frequency and decay 
rate /^m- We allow x to couple to the two optical modes 
through the coupling constants ^f,s; we shall make the 
simplifying assumption that, since ^f,s is proportional to 
the respective optical frequency, = 2^f- Thus, we can 
write the Hamiltonian ^ as a sum of four terms 



H — i^free + ^diss + ^OM + ^SHG 



(2) 



with the free Hamiltonian (we take units such that h = 1 
throughout the paper) 



^free= J duj uj ala^^ ^ 



liijUJca]dj^—-{x +j9 ), (3) 

J=F,S 

where /i^ = 2 and /i^? = 1, the dissipation Hamiltonian 



M.diss 



, (4) 



where we leave the mechanical dissipation Hamiltonian 
undefined, the optomechanical Hamiltonian Hqm = 
— J2j=F s 9j^]^j^ and, finally, the SHG Hamiltonian [10] 



SHG 



iX {alfas - (^f)^^s • 



(5) 



In the expression for i/^diss we defined two frequency 
ranges flF,s^ which define the bath modes through which 
the two modes are damped. Given the very large separa- 
tion in frequency between the two optical modes, we can 
justify considering and Qs as non-overlapping with- 
out violating the requirement that these two frequency 
ranges must be very large compared to hz^^s that is nec- 
essary to ensure Markovian dynamics. If the two fre- 
quency ranges were to have a significant overlap, e.g., in 
the case of two optical modes spaced by a mechanical fre- 
quency that exist in different resonators but are coupled 
to the same bath, then one must be careful to use cor- 
rectly modified input-output relations and equations of 
motion. In our case, the standard input-output relations 



will be held valid. In the following we will obtain the 
linearized equations of motion for the system described 
by this Hamiltonian. 

As a first step we derive the Heisenberg equations of 
motion for as,F^ which read: 



as 



(zA — a>:f)^f ^ 
{2iA - Ks)cts 



- x(^f)^ - 




where ftp^g) is the input field coupled to the funda- 
mental mode (the second harmonic) having zero mean 
and two-time correlation function (%(g) (^)^p(g) (^0) = 
S{t — t') and A = cjp — is the detuning of the 
driving field from cavity resonance. Eqs. (6) are writ- 
ten in a frame rotating with the optical modes, i.e., in 
an interaction picture with respect to the Hamiltonian 
X]j=s F (/q ^lo^oj ^]%) • Similarly, for the me- 

chanical mode we have x = uj^p and 



P = -UJmX - 2f^rnP - V2^m^ + ^F^p^F + ^S^^S, (7) 

where we used a Brownian-motion-type damping 
model [11]. The self-adjoint Langevin force ^ is zero- 
mean and (assuming the high-temperature limit) delta- 
correlated as {i{t)i{t')) = {2nth^l)S{t-t'). We now lin- 
earize the equations of motion by considering a pumping 
field of large intensity. Under these conditions, both the 
field modes of the toroid would be macroscopically pop- 
ulated. We are then allowed to take aF,s = ^f,^ + ^^^f,S', 
where dF,s = {^f,s) is the (large) mean amplitude of 
each operator and Sd is its fiuctuation around such aver- 
age. For simplicity of notation, we shall drop the 'J' in 
the operator fiuct nations. We thus write (j = F, S) 



[i(/ijA + g^x) — /^jjap + ig^a^x - 



(8) 



with Of = 2x{apds + agap) and Og = — 2xttF^F- De- 
spite the relation linking g^ and ^g, in what follows we 
shall continue to use both symbols for clarity. For the 
mechanical modes, the equation of motion for x remains 
unchanged, while that for momentum becomes 



P 



aX - 2k^P - V^Kmi - 



■E 

J=F,S 



Assuming, for now, the existence of a steady state we 
find the following equation relating the amplitude of the 
fundamental mode and that of the input noise 



(10) 



where we have taken (j) is the phase of ap. The input 
field is assumed to be in a monochromatic coherent state 
characterized by the (real) amplitude ain = ap . For the 
parameters used throughout this paper, the first term on 
the left-hand side in the above equation can be safely 
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(a) Reductions, Qm = 5970 (b) Bipartitions, Qm = 5970 (c) Reductions, Qm = 597000 (d) Bipartitions, Qm = 597000 

FIG. 2. (Color online) Regions of (a) 1-1-mode ('reduction') entanglement, or (b) of 1-2-mode ('bipartition') entanglement. In 
subfigure (a) the dashed curve encompasses the (red) region where the two-mode reduction consisting of the second-harmonic 
optical mode (system S) and the mechanical mode M is not separable, the dotted curve encompasses the (green) region where 
the fundamental optical mode (system F) and the mechanical mode are not separable, and the dashed-dotted curve delineates 
the (blue) region where the two optical modes are entangled. A similar explanation holds for subfigure (b), but with the 
dashed curve corresponding to entanglement between the fundamental optical mode and the 2-mode system formed by the 
second-harmonic mode and the mechanical mode. The dotted curve delineates inseparability of the second-harmonic from the 
2-mode system formed by the other two modes, and likewise the dashed-dotted curve bounds the region where the mechanical 
mode is entangled with the 2-mode system consisting of the optical modes. The solid curve encloses the region of instability, 
(c) and (d) are similar to (a) and (b), respectively, but with a larger mechanical Q-factor. Because of the smaller mechanical 
decay rate, the region of instability starts at significantly lower powers, especially on the blue-detuned side of the figure. 



neglected, whereupon the equation can easily be solved 
to obtain ap. For the second harmonic, we get 



as 



2i</) 



2iA - K.S 



(11) 



while X = p = 0. The linearization of the equations of 
motion makes the dynamics Gaussian so that any initial 
Gaussian will remain such at any instant of time [12]. We 
now introduce the quadrature vector 



(12) 



where we have defined the optical quadrature operators 
f j = -l=(ai + al) and = -^{^] ~ ^i) mode j = F, S. 
Similar definitions hold for the input-field operators. It is 
worth bearing in mind that all our operators are functions 
of time and we have simply dropped the label t for 
conciseness of notation. With these definitions, the first 
moment of R is zero and any Gaussian state of the system 
is thus fully characterized by the covariance matrix V = 
( {R^ R)-\-{R<^R)^)/2. The equations of motion derived 

above can be concisely written diS R = A - R-\- i^in, with 
the input noise vector 



R\r} 



(13) 

The drift matrix A can be explicitly determined and de- 
pends on the set of parameters characterizing the dy- 
namics of the three- mode system addressed here. Its ex- 
pression is too lengthy to be reported here and is thus 
deferred to the Appendix. A close inspection of the form 
of A reveals that, by assuming ain G M and introducing 
the rescaled parameters a = ^f/(a/2x) = X^f, 

the drift matrix is a universal expression of Qf/x- Thus, 
for a fixed value of /3 the nature of the dynamics is de- 
termined solely by the ratio of the coupling constants. 



For small a, the interaction is dominated by the SHG 
process. Conversely, for a 1, the dynamics resembles 
closely that of a standard optomechanics model with two 
fields [13]. The dynamical equations should be stable in 
order for a steady state to exist. This is assured if the real 
part of the spectrum of A is negative, in which case the 
system will tend to a stationary state characterized by 
the covariance matrix that solves the Lyapunov equation 
A • V + V • A^ + D = with the input-noise matrix 

D 5{t-t')=\ [{RUt)®RUt'))+{RUt)®Rin{t')f] . 

(14) 

The stability condition for A can be rephrased more 
formally in terms of the Routh-Hurwitz criterion [14], 
which we have used in our quantitative characterization 
of the dynamics. We note in passing that the Lya- 
punov equation above has a concise analytic solution 
for V, as reported in Ref. [14]. The covariance ma- 
trix encompasses the full information on the system at 
hand. Here, we shall be interested in the entanglement- 
sharing properties of the three modes. In order to 
demonstrate the occurrence of genuine multipartite en- 
tanglement, we rely on the criterion based on nega- 
tivity of partial transpose (NPT) [15-17] and we will 
make use of the logarithmic negativity as an entangle- 
ment quantifier [18]. For a bipartition consisting of sub- 
systems A and B (A, B = F,S,M), this is defined as 



E 



A|B 



max[0, — ln(2 z>_^/e)] , where {z^-,/e} is the 
set of symplectic eigenvalues of the covariance matrix 
associated with the partially transposed states of the 
system such that |z^-,/c| < 1/2 [19]. If either A or B 
are single- mode subsystems, k = 1 regardless of the 
number of modes comprised in B or A. Moreover, al- 
though this is not the case in general, in this situation 
the NPT criterion is a necessary and sufficient condi- 
tion for inseparability of pure and mixed Gaussian states 
alike. In what follows, we characterize the entanglement 
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FIG. 3. (Color online) Tripartite logarithmic negativity, E)^^ ^ for five different input powers and both low- and high-Qm cases. 



structure in both reduced two-mode states and bipar- 
tite one- versus- two-mode ones. In order to do that, we 
use numerical values for the various constants entering 
the model that reflect the state-of-the-art of recent ex- 
periments. The fundamental wavelength is chosen to be 
1554 nm [8], at which an input power Pin = 1 M-W cor- 
responds to 10^ s ^. Moreover, we set cjm = 
27rx70MHz,/^m = 27rx5.9kHz (Q^ 5970), /^f = /^s = 
27r X 7MHz,^/F = Qs/^ = 27r x 1.2kHz,Tenv = 0.8 K [20], 
and X = '^OO Hz, which is within a factor of 2 of what has 
been observed in Ref. [8]. 

For these choices. Fig. 2 summarizes both the entan- 
glement in one-versus-one-mode reduced states and one- 
versus-two-mode situations. Entanglement is analyzed in 
the A-Pin parameter space. Pin being the input power. 
While all-optical entanglement (i.e., the entanglement 
within the reduction involving only the S and F sub- 
systems) exists in a narrow strip around A c:^ and is 
maximum on resonance, the mechanical mode is entan- 
gled more strongly with S near A = ± cjm, and with F 
close to A = ± At low (yet still quite sizeable) 

values of the mechanical quality factor, the region corre- 
sponding to A < is largely associated with separability 
of any two-mode reduction, except the narrow strip at 
A :^ mentioned above, which witnesses the fact that, 
in these conditions, the direct nonlinear coupling between 
the optical modes overcomes any entangling power of 



TABLE L Calculated logarithmic negativities for A = — cjm 
and Pin = 0.27W in Figs. 2(c) and (d). The rest of the 
parameters as in the body of the paper, with Qm = 597000. 



Reduction a\b 




Bipartition a\bc 


jpa 1 he 


S|M 


0.10 


F|SM 


0.44 


F M 


0.42 


S FM 


0.15 


FS 


0.01 


M FS 


0.45 



the optomechanical mechanism. Larger values of Qm, on 
the other hand, give rise to non-negligible areas of (even 
strong) optomechanical entanglement involving both the 
M-S pair and the M-F one. This is clearly shown in 
Fig. 2(c) where, remarkably, we find that even the all- 
optical entanglement is affected, spreading quite consid- 
erably in regions where, at lower Qm, we had E'ff = 0. 
This result might be interpreted as arising from the indi- 
rect coupling between the two optical modes, whose inter- 
action is ruled not only by their direct nonlinear coupling 
but also by a detuning-dependent effective one mediated 
(quasi-coherently) by the mechanical mode. 

Also quite interesting is the behavior of the one-vs- 
two-mode entanglement. An investigation on these con- 
figurations is relevant in order to characterize the multi- 
partite entanglement being possibly shared by the three 
subsystems. Indeed, based on the classification provided 
by Giedke et al. [21], the simultaneous inseparability of 
the three possible one-versus-two-mode bipartitions in 
a three-mode system implies the existence of genuine 
tripartite entanglement. Likewise, the state is /c-mode 
biseparable if there are k one-versus-two-mode biparti- 
tions with respect to which the state of the system is 
separable. Fig. 2(b) and Fig. 2(d) show the rich structure 
of entanglement sharing that is exhibited by our model. 

One-versus-two-mode entanglement turns out to be, 
in general, much more robust (and larger) with respect 
to noise affecting the system than the entanglement in 
any two-mode reduction, a feature that has already been 
shown in other optomechanics-related investigations [22] . 
In Fig. 2(b), regions of full three- mode inseparability are 
shown even for a relatively low-quality mechanical os- 
cillator. However, in this case, noise affecting the sys- 
tem through the mechanics is too strong to allow for 
much overlap between regions of three-mode insepara- 
bility. Indeed, by increasing the mechanical quality fac- 
tor by a factor of 100, giving the results in Fig. 2(c), 
the overlap between regions of one-versus-two-mode en- 
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FIG. 4. (Color online) Similar to Fig. 2, but varying the second- harmonic generation rate on the vertical axis, (a) and (b) have 
= 5970, whereas (c) and (d) have = 597000. (xo = 27r x 700 Hz, Pin = 10"^ W.) 



tanglement increases significantly, covering virtually the 
whole stability area shown in the figure. In passing, we 
mention that we have applied a multipartite entangle- 
ment witness for continuous-variable states, MultiWit, 
that was developed in Ref. [23] using semi-definite opti- 
mization methods. The use of this instrument has con- 
firmed the genuinely tripartite nature of the entangle- 
ment at hand in the regions of overlap among the three 
regions of one-versus-two-mode inseparability, which ex- 
cludes the possibility of having generalized three-mode 
biseparable states. As a quantitative illustration for the 
high-Qm case, in Table I we give the entanglement in 
any reduction and bipartition that can be singled-out in 
our problem, taking the values of the parameters listed 
above and choosing A = —oOrn with Pin = 0.27 W. In or- 
der to complete our assessment, we have determined the 
degree of genuine tripartite entanglement across interest- 
ing regions in the full-inseparability areas. As a quanti- 
tative estimator, we have used the tripartite logarithmic 
negativity Ej^\ which is a proper entanglement mono- 
tone [24] . The results of this study are shown in Fig. 3 for 
both the low- and high-Qm cases and increasing optical 
input powers. The qualitative differences in the behavior 
of the tripartite entanglement is very marked at large in- 
put power: whilst at low mechanical quality factors the 
high-power tripartite entanglement is null, it extends for 
most of the region A G [— cJm, 0] at high mechanical qual- 
ity factor, therefore leaving us with much room for ma- 
neuvering in the space of entangled three-mode states. 
Let us finally explore the competition between the SHG 
process and the optomechanics in our model. As x grows 
in Fig. 4(a) and Fig. 4(c), the entanglement for the F|M 
and S|M reductions decreases, whereas that for F|S cov- 
ers an ever-larger area of the parameter space. The be- 
havior of the entanglement in the bipartitions. Fig. 4(b) 
and Fig. 4(d), is similar and also easily understood on 
an intuitive basis: a bigger x leads to larger regions of 
entanglement for the two bipartitions that involve one of 
the optical mode on its own (i.e., F|MS and S|MF), but a 
contraction in the parameter space where entanglement 
in the M|FS bipartition is observed. 



II. INFERRING THE STATE OF THE SYSTEM 

The inference of the full state of an optomechanical sys- 
tem is a major practical challenge [13], mainly due to the 
fact that the mechanical quadratures are not directly ac- 
cessible to an experiment. Here we propose a technique, 
which requires the use of the system drawn schematically 
in Fig. 5, that allows us to infer the mechanical quadra- 
tures indirectly. Initially, we assume that the homodyne 
detectors needed in the scheme have infinite bandwidth; 
we shall account for the finite bandwidth of any realistic 
apparatus later on. Given such a system, one has ac- 
cess to the four input (two each for the fundamental and 
the second harmonic) and four output quadratures. One 
can then infer the intra-cavity optical modes using the 
input-output relation 



(15) 



and similarly for the rest of the quadratures. It is clear 
that very good characterization of the system, includ- 
ing knowledge of all the coupling constants and the ef- 
fect of vacuum input noise [13] is required to infer the 
intra-cavity quadratures accurately. Once the four opti- 
cal intra-cavity quadratures are known, it is natural to 
ask how the mechanical quadratures can be inferred. In- 
deed, the key to our proposal is noticing that the inferred 
quadratures are obtained as a time trace. One can there- 
fore make direct use of the equations of motion to obtain 
a time trace for the mechanical quadrature x and p. The 
covariances of these inferred quadratures can finally be 
used to build an inferred covariance matrix, V. In the 
limit of infinite detector bandwidth, V = V. 

We use a simple model for including the effects of a fi- 
nite detector bandwidth. Let the bandwidth of the detec- 
tor be r. Then the point-spread function of the detector 
is taken to be 



fity- 



e{t) - e{t 





1 

T 

J_ 

2r 



for t <0 OT t> r 

for t e (0,r) 

for t = or t = r 



(16) 

where B(t) is the Heaviside step function. The label t is 
understood as the time at which the measurement was 
performed. Normalization requires that f{t)dt = 1. 



6 



Output fields 



^out 
I ^out 



^out 



Input fields 



Ps 



FIG. 5. Schematic of detection system. Part of the input field 
is used as the local oscillator of two homodyne detectors (this 
necessitates adding the second harmonic onto the local oscil- 
lator), one recording the field quadratures before the toroid, 
and one after. 



For each operator a{t) we assign an inferred operator 
d{t) = {f^d){t) = f{t — s)d{s) ds as the convolution of 
f{t) with d{t). It can easily be shown that ^(/*^)(^) = 
(/ * d){t). This means that we can infer the value of 
any d{t) by calculating the time derivative of the inferred 
~a{t). 

With this at hand, we can finally show that the inferred 
covariance matrix, at the steady-state, is given by 



V= (AT)"'{(e^^-l)V(e^--l) 



Jo 



1]D[, 



Mrs) 



ifds 



}( 



At) 



(17) 

Modern homodyne detectors can operate with a band- 
width of the order of 10 GHz, which is much larger than 
typical values of ujm in the micro-mechanical domain. 
Therefore, it is understood that r is by far the short- 
est timescale of the system. It then suffices to expand V 
to first order in r, which gives 



V = V- 



(18) 



These expressions hint at the tantalizing possibility that 
by increasing r electronically one could deduce the value 
of V for vanishing r, and therefore infer V itself. In our 
numerical exploration, the fidelity [25] for the inference of 
the mechanical mode using this method was above 99% 
when r corresponded to a bandwidth of 500 MHz. 



III. CONCLUSIONS AND OUTLOOK 

We have presented a system that combines a nonlin- 
ear optical process with optomechanics in a very natural 
manner. Its monolithic design makes it very attractive 
for experimental, or even technological, applications. In- 
deed, the system we presented is based on technology 
that is inherently compatible with integration on opto- 
electronic chips. Our investigation concentrated on the 
dynamics of the system, but we also addressed the prob- 
lem of the actual detection of the intra-cavity state by 
outlining a method involving homodyning all the input 
and output quadratures to infer the covariance matrix 
for the three intra-cavity modes. Lastly, a numerical ex- 
ample using constants from recent experiments was used 
to illustrate the feasibility of observing these effects in a 
realistic system. 

Looking further ahead, one can envisage several of these 
structures sharing a common photonic "bus" whose func- 
tion is to populate the optical mode of each toroid at the 
fundamental frequency. The SHG process and optome- 
chanics could then be used to create an entangled state 
of the mechanics with the second-harmonic field, thereby 
generating for each structure an optical field, which can 
be routed out of the photonic bus without losses, due 
to the large separation in frequencies, that is entangled 
to the mechanical mode of the toroid. We remark that 
this system lends itself naturally to the distribution and 
certification of optomechanical entanglement as per the 
protocol recently proposed in Ref. [26]. 
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APPENDIX 

Here we provide the explicit form of the drift matrix 
A for our problem, which reads 





-A + 2x4 


2xa^F 


2xa^ 


-v^^f4 





A + 2x4 


-f^F - 2X^8 


-2x4 


2x^F 


V^^F^F 







2x4 


-i^s 


-2A 







-2x^F 


-2x^F 


2A 


-^s 


V^gsal 






















V2^F^F 




V^gs^s 




-UJm 


-2n, 



(19) 
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with aj = aj + m- {j = F, S). As done in the main body of the paper, we assume that ap is real and ^fs = 2^f- This 
can always be done by dropping the corresponding assumption on ain and choosing its phase appropriately. We then 
define the parameters ol — Qy j (a/2 x) /3 = x^f rewrite the drift matrix as 





-A- 


4A 


4A (Yl 


4A2+m:| 
2m:s 




4A2+/.C2 


-2/3 












-2/3 










2a/3 








2/3 


2A 


4A2 + h:2 



a/32 





2/3 
-2A 





8A 
4A2+/^2 



a/32 





2a/3 



8A 



4A2+4«/^' 
"4A2+4 










— 2/^ri 



showing that the steady-state of the system is a universal /3 the dynamics is determined solely by the value of the 
function of a and /3. In particular, at a fixed value for ratio of the coupling constants a. 
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